function out = int_r( k, r, t, s_, t_, r_)
%INT_R Summary of this function goes here
%   Detailed explanation goes here

global h_theta;

out = r_.^2 .* trpzInt(@(theta)int_th(k, r, t, s_, t_, r_, theta), 0, pi, h_theta);

end

